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Abstract. In 1991, David Gale and Raphael Robinson, building on explorations carried out 
by Michael Somos in the 1980s, introduced a three-parameter family of rational recurrence 
relations, each of which (with suitable initial conditions) appeared to give rise to a sequence 
of integers, even though a priori the recurrence might produce non-integral rational numbers. 
Throughout the '90s, proofs of integrality were known only for individual special cases. In 
the early '00s, Sergey Fomin and Andrei Zelevinsky proved Gale and Robinson's integrality 
conjecture. They actually proved much more, and in particular, that certain bivariate ratio- 
nal functions that generalize Gale-Robinson numbers are actually polynomials with integer 
coefficients. However, their proof did not offer any enumerative interpretation of the Gale- 
Robinson numbers/polynomials. Here we provide such an interpretation in the setting of 
perfect matchings of graphs, which makes integrality/polynomiality obvious. Moreover, this 
interpretation implies that the coefficients of the Gale-Robinson polynomials are positive, as 
Fomin and Zelevinsky conjectured. 



In memory of David Gale, 1921-2008 
1. Introduction 

Linear recurrences are ubiquitous in combinatorics, as part of a broad general framework that 
is well-studied and well-understood; in particular, many combinatorially-defmed sequences can 
be seen on general principles to satisfy linear recurrences (see [26]), and conversely, when an 
integer sequence is known to satisfy a linear recurrence it is often possible to reverse-engineer 
a combinatorial interpretation for the sequence (see [I] and references therein for a general 
discussion, and [U Chapter 3] for specific examples). In contrast, rational recurrences such as 

s{n) = (s(n - l)s(n - 3) + s(n - 2) 2 )/s(n - 4), 

which we prefer to write in the form 

s(n)s(n - 4) = s(n - l)s(n - 3) + s(n - 2) 2 , 

are encountered far less often, and there is no simple general theory that describes the solutions to 
such recurrences or relates those solutions to combinatorial structures. The particular rational 
recurrence relation given above is the Somos-4 recurrence, and is part of a general family of 
recurrences introduced by Michael Somos: 

s(n)s(n -k) = s(n - l)s(n - k + 1) + s(n - 2)s(n - k + 2) + h s(n - [k/2\)s(n - \k/2\). 

If one puts s(0) = s(l) = • • • = s(k — 1) = 1 and defines subsequent terms using the Somos- 
k recurrence, then one gets a sequence of rational numbers which for the values k = 4, 5, 6, 7 
is actually a sequence of integers. (Sequences Somos-4 through Somos-7 are entries A006720 
through A006723 in [24].) Although integer sequences satisfying such recurrences have received 
a fair bit of attention in the past few years, until recently algebra remained one step ahead of 
combinatorics, and there was no enumerative interpretation of these integer sequences. (For 



links related to Somos sequences, see http://jamespropp.org/somos.html ) 
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Inspired by the work of Somos, David Gale and Raphael Robinson [HI Qj] considered se- 
quences given by recurrences of the form 

a(n)a(n — m) = a(n — i)a(n — j) + a(n — k)a(n — £), 

with initial conditions a(0) = a(l) = • • ■ = aim — 1) = 1, where m = i + j = k + 1 . We call 
this the three-term Gale-Robinson recurrence^. The Somos-4 and Somos-5 recurrences are the 
special cases where (i, j, k, £) is equal to (3, 1,2,2) and (4, 1, 3, 2) respectively. Gale and Robinson 
conjectured that for all integers i, j, k,£ > with i + j = k + £ = m, the sequence a(0), a(l), . . . 
determined by this recurrence has all its terms given by integers. About ten years later, this was 
proved algebraically in an influential paper by Fomin and Zelevinsky 

1.1. Contents 

In this paper, we first give a combinatorial proof of the integrality of the three-term Gale- 
Robinson sequences. The integrality comes as a side-effect of producing a combinatorial in- 
terpretation of those sequences. Specifically, we construct a sequence of graphs P{n]i,j,k,£) 
(n > 0) and prove in Theorem [9] that the nth graph in the sequence has a{n) (perfect) match- 
ings. Our graphs, which we call pinecones, generalize the well-known Aztec diamond graphs, 
which are the matchings graphs for the Gale-Robinson sequence 1, 1, 2, 8, 64, 1024, ... in which 
i = j = k = £ = 1. A more generic example of a pinecone is shown in Figure [TJ All pinecones 
are subgraphs of the square grid. 



Figure 1. The pinecone P(25; 6, 2, 5, 3). Its matching number is a(25), where 
a(n) is the Gale-Robinson sequence associated with (i,j,k,£) = (6,2,5,3). 

We give two ways to construct pinecones for the Gale- Robinson sequences: a recursive method 
(see Figure [TT] and the surrounding text) that constructs the graph P(n; k, £) in terms of the 
smaller graphs P (n'\ i,j, k, £) with nl < n, and a direct method (see Formula ([2]) in Section[3]) that 
allows one to construct the graph P(n; k, £) immediately. The heart of our proof is the demon- 
stration that if one defines a(n) as the number of perfect matchings of P(n) = P(n;i,j,k,£), 
the sequence a(0), a(l), a(2), ... satisfies the Gale-Robinson recurrence. This fact, in combina- 
tion with a simple check that a(0) = a(l) = • • • = a(m — 1) = 1, gives an immediate inductive 
validation of our claim that P(n) has a(n) perfect matchings for all n, which yields additionally 
the integrality of a(n). 

General pinecones are defined in Section [2j where we also explain how to compute induc- 
tively their matching number via Kuo's condensation lemma [T7]. In Section O we describe 
how to associate a sequence of pinecones to a Gale-Robinson sequence, and observe that for 
these pinecones, the condensation lemma specializes precisely to the Gale-Robinson recurrence. 
Indeed, the recursive method of constructing pinecones, in combination with Kuo's condensa- 
tion lemma, gives combinatorial meaning to the different terms a{n\)a{n2) of the Gale-Robinson 
recurrence. 

In Section 21 we refine our argument to prove that the sequence p(n) = p(n; u>, z) defined by 

p(n)p(n — m) = w p(n — i)p(n — j) + zp(n — k)p(n — £), 

^Gale and Robinson also considered recurrences of the form a(n)a(n — m) = a{n — g)a(n — h) + a{n — i)a{n — 
j) + a(n — k)a{n — £) for suitable values of g, h,i, j, k, £,m, but such four-term Gale- Robinson recurrences will not 
be our main concern here. 
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with i+j = k + i = m and p(0) = p(l) = ■ ■ ■ = p(m — 1) = 1, is a sequence of polynomials 
in w and z with nonnegative integer coefficients. More precisely, we prove in Theorem l20l that 
p(n;u 2 ,v 2 ) counts perfect matchings of the pinecone P(n;i,j,k,£) by the number of special 
horizontal edges (the exponent of the variable u) and the number of vertical edges (the exponent 
of the variable v). The fact that p(n) is a polynomial with coefficients in Z was proved in [TT] . 
but no combinatorial explanation was given and the non-negativity of the coefficients was left 
open. 

1.2. Strategy, and connections with previous work 

For much of the work in this paper, we share precedence with the students in the NSF- 
funded program REACH (Research Experiences in Algebraic Combinatorics at Harvard), led 
by James Propp, whose permanent archive is on the web at http://jamespropp.org/reach/, 
A paper by one of these students, David Speyer [25], introduced a very flexible framework 
(the "crosses and wrenches method") that, starting from a recurrence relation of a certain type, 
constructs a sequence of graphs whose matching numbers satisfy the given recurrence. This 
framework includes the three-term Gale- Robinson recurrences, and thus yields a combinatorial 
proof of the integrality of the associated sequences. This extends to a proof that the bivariate 
Gale- Robinson polynomials mentioned above are indeed polynomials, and have non-negative 
coefficients. One difference with our paper is that Speyer's graphs are only described explicitly 
for Somos-4 and Somos-5 sequences, whereas our construction is explicit for any Gale-Robinson 
sequence. Moreover, the description of our graphs as subgraphs of the square grid looks more 
regular, and may be useful to study limit shapes of random perfect matchings such as the perfect 
matchings shown in Figure [19l 

Let us mention that shortly after Speyer did his work on perfect matchings, he and his fellow 
REACH-participant Gabriel Carroll did for four-term Gale-Robinson recurrences what Speyer 
had done for three-term Gale-Robinson recurrences, by introducing new objects called "groves" 
to take the place of perfect matchings [6]. Carroll and Speyer's work gives, as two special cases, 
combinatorial proofs of the integrality of Somos-6 and Somos-7. 

The strategies that led to Speyer's article [25] and to the present article are not entirely 
independent; each made use of Propp 's prior construction of a suitable perturbed Gale- Robinson 
recurrence, which we explain next. The explanation will mostly be of interest to researchers 
seeking to apply similar techniques to other problems; others may want to skip the rest of the 
introduction. 

Suppose we perturb a three-term Gale-Robinson recurrence by replacing the singly-indexed 
Gale- Robinson number a(n) by a triply-indexed quantity A(n,p, q) satisfying the perturbed 
recurrence 

A(n,p 7 q)A(n - m 7 p, q) = A(n — i,p—l, q)A(n — j,p + 1, q) + A(n - k,p,q + l)A(n - £,p, q— 1). 

(This choice of perturbation is not as special as it looks: all that matters is that the pairs 
(—1,0), (1,0), (0, 1), (0, —1) that describe the perturbations of the second and third coordinates 
in the four index-triples on the right-hand side, viewed as points in the plane, form a non- 
degenerate centrally-symmetric parallelogram. Choosing a different centrally-symmetric paral- 
lelogram is tantamount to a simple re-indexing of the recurrence.) If we take as our initial 
conditions A(n,p, q) = x n ,p, q for all n between and m — 1 and p, q arbitrary, with (formal) 
indeterminates x ntP . q , then each A(n,p, q) with n > m can be expressed as a rational function of 
these indeterminates. It should be emphasized here that for all n,p, q, r, s, the rational functions 
A(n,p, q) and A(n, r, s) are the same function up to re-indexing of the indeterminates. 

Propp conjectured that each A(n,p, q) is a Laurent polynomial in some finite subset of the 
(infinitely many) indeterminates x n ^ s , with integer coefficients; that is, each A(n,p, q) is an 
element of Z[a^* 8 ]. This was subsequently proved by Fomin and Zelevinsky [11]. Note that if 
one sets all the indeterminates x n , r ,s equal to 1, the Laurent polynomials A(n,p, q) specialize 
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to the Gale- Robinson numbers a(n). Propp conjectured that each coefficient in each such Lau- 
rent polynomial is positive (a fact that is not proved by Fomin and Zelevinsky's method) and 
furthermore is equal to 1. 

Propp knew that in the case i = j = k = £=l, the Laurent polynomials A(n,p,q) can be 
interpreted as multivariate matching polynomials of suitable graphs, namely, the Aztec diamond 
graphs. (See Subsection 12.11 for a definition of matching polynomials.) Indeed, David Robbins 
had studied the three-parameter "perturbed recurrence" in this case, on account of its relation 
to the study of determinants, and had shown (with Rumsey) [22j that the associated rational 
functions are Laurent polynomials. (For more background on this connection with determinants, 
see [5].) The work by Elides, Kuperberg, Larsen, and Propp [10] had shown that the monomials 
in these Laurent polynomials correspond to perfect matchings of Aztec diamond graphs. So it 
was natural to hope that this correspondence could be extended to the Gale-Robinson family of 
recurrences. 

It should be acknowledged here that the idea behind the specific triply-indexed perturba- 
tion A(n,p, q) of the Gale- Robinson sequence that proved so fruitful came from an article of 
Zabrodin [28] that was brought to Propp's attention by Rick Kenyon. This article led Propp 
to think that the recurrence studied by Robbins should be considered a special case of the "dis- 
crete bilinear Hirota equation", or "octahedron equation", and that other recurrences such as 
the Gale-Robinson recurrence should likewise be considered in the context of the octahedron 
equation. 

What the REACH students were able to do, after diligent examination of the Laurent poly- 
nomials A(n,p,q), is view those Laurent polynomials as multivariate matching polynomials of 
suitable graphs. Bousquet-Melou and West, independently, did the same for small values of n, 
until they were able to extrapolate these examples to the generic form of the graphs, which 
became the pinecones of this paper. 

There is a general strategy here for reverse-engineering combinatorial interpretations of 
algebraically-defined sequences of numbers: add sufficiently many extra variables so that the 
numbers become Laurent polynomials in which every coefficient equals 1. For another applica- 
tion of this reverse-engineering method (in the context of Markoff numbers and frieze patterns) , 
see Q3]. 

2. Perfect matchings of pinecones 

In this section we define a family of subgraphs of the square lattice, which we call pinecones. 
Then we prove that the number of perfect matchings of a pinecone can be computed inductively 
in terms of the number of perfect matchings of five of its sub-pinecones. 

2.1. Preliminaries 

To begin with, let us recall some terminology about graphs. A (simple) graph G is an ordered 
pair (V, E) where V is a finite set of vertices, and E, the set of edges, is a collection of 2-element 
subsets of V. The degree of a vertex v is the number of edges in E containing v. A subgraph of G 
is a graph H = (V, E') such that V C V and E' C E. If, in addition, V = V, we say that H is a 
spanning subgraph of G. The intersection of two graphs G = (V, E) and H = (V , E') is the graph 
GnH = (Vr\V',Er\E'), and the union of the two graphs is the graph GUH = (Vl)V',El)E'). 
Given two graphs G = (V, E) and H = (V, E'), we denote by G \ if the subgraph (V",E"), 
where V" = V \ V and E" is the set of edges of E \ E' having both endpoints in V" . 

A perfect matching of a graph G = (V, E) is a subset E' of E such that every vertex of V 
belongs to exactly one edge of E'. We will sometimes omit the word "perfect" and refer to perfect 
matchings as simply "matchings". The matching number of G, denoted by m(G), is the number 
of perfect matchings of G. More generally, we shall often consider the set E of edges as a set of 
commuting indeterminates, and associate with a (perfect) matching E' the product of the edges 
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it contains. The matching polynomial of G is thus defined to be 

M (G) :=E II e ' 

E' e£E' 

where the sum runs over all perfect matchings E' of G. If we replace every e that occurs in 
this sum-of-products by a non-negative integer n e , then this expression becomes a non-negative 
integer, namely, the number of perfect matchings of the multigraph in which there are n e edges 
joining the vertices x and y for all e = {x, y} in E (and no edges joining x and y if {x, y} is not 
in E). In particular, if each n e is set equal to or 1, then the matching polynomial becomes 
the number of perfect matchings of the subgraph of G consisting of precisely those edges e for 
which n e = 1. 

9 



n 




s 



Figure 2. An Aztec diamond graph of width 9, and one of its perfect matchings. 

2.2. Aztec diamonds graphs 

The pinecones considered in this paper are certain subgraphs of the square lattice. The 
most regular of them are the (Aztec) diamond graphs, which are the duals of the so-called 
Aztec diamonds, which were first studied in detail in [10]. A diamond graph of width 2k — 1 is 
obtained by taking consecutive rows of squares, of length 1, 3, ... , 2k— 3, 2k— 1, 2k— 3, . . . , 3, 1 and 
stacking them from top to bottom, with the middle squares in all the rows lining up vertically, 
as illustrated by Figured! Let A be a diamond graph of width 2k — 1. Let An be the diamond 
graph of width 2k — 3 obtained by deleting the leftmost and rightmost squares of A as well as 
the two lowest squares of each of the remaining 2k — 3 columns of A. We call An the North 
sub-diamond of A. Define similarly the South, West and East sub-diamonds of A, denoted by 
As, Aw and Ae- Finally, let Ac be the central sub-diamond of A of width 2k — 5 (Figure[3]). The 
following result is a reformulation of Kuo's condensation theorem for Aztec diamond graphs [T7] . 

Theorem 1 (Condensation for diamonds graphs). The matching polynomial of a diamond 
graph A is related to the matching polynomials of its sub- diamonds by 

M(A)M(A C ) = nsM(A w )M(A E ) + ewM{A N )M(A s ), 

where n, s,w, and e denote respectively the top (resp. bottom, westmost, eastmost) edge of A (see 
Figure^). 

In particular, if a(n) (with n > 2) denotes the matching number of a diamond graph of width 
2n — 3, then 

a(n)a(n - 2) = 2a(n - l) 2 
for all n > 2, provided we adopt the initial conditions a(0) = a(l) = 1. This shows that a(n) 
is the three-term Gale- Robinson sequence associated with i — j = k = £=l, and implies 
a(n) = 2ft). 

The condensation theorem we shall prove for pinecones appears as a generalization of the 
condensation theorem for diamond graphs. But it can actually also be seen as a specialization 
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Figure 3. The five sub-diamonds of a diamond graph of width 9. 



of it, and this is the point of view we adopt in this paper. The key idea is to forbid certain edges 
in the matchings. 

Corollary 2. Let A be a diamond graph, and let G be a spanning subgraph of A, containing the 
edges n, s, w and e. Let Gn — G (~l An , and define Gs, Gw, Ge and Gc similarly. Then 

M{G)M(G C ) = nsM{Gw)M(G E ) + ewM(G N )M(G s ). 

Proof. Since G is a spanning subgraph of A, every perfect matching of G is a perfect matching 
of A. Hence the matching polynomial M(G) is simply obtained by setting a = in M(A), for 
every edge a that belongs to A but not to G. The same property relates M(Gjv) and M(An), 
and so on. Consequently, Corollary [2] is simply obtained by setting a = in Theorem [TJ for 
every edge a that belongs to A but not to G. □ 

2.3. PlNECONES: DEFINITIONS 

A standard pinecone of width 2k — 1 is a subgraph P = (V, E) of the square lattice satisfying 
the three following conditions, illustrated by Figure Hla: 

1. The horizontal edges form i + j + 1 segments of odd length, starting from the points 
(0, 1), (1, 2) . . . , (i - 1, i) and (0, 0), (1, -1), . . . , for some i > 1, j > 0. Moreover, 
if L m denotes the length of the segment lying at ordinate m, then 

L_j < • • • < £_i < L Q = 2k - 1 = ii > L 2 > ■ ■ ■ > 

2. The set of vertices V is the set of vertices of the square lattice that are incident to the 
above horizontal edges. 
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3. Let e = {(a, b), (a, b + 1)} be a vertical edge of the square lattice joining two vertices of 
V. If a + b is even, then e belongs to the set of edges E, and we say that e is an even 
edge of P. Otherwise, e may belong to E, or not (Figure HJa), and we call e a (present 
or absent) odd edge. 

The leftmost vertices of a standard pinecone are always (0, 0) and (0, 1). However, sometimes it 
is convenient to consider graphs obtained by shifting such a graph to a different location in the 
two-dimensional lattice. We will call such a graph a transplanted pinecone. In a transplanted 
pinecone, the leftmost vertices are (a, b) and (a, 6+1), where a + b is even. In some cases, where 
the distinction between standard and transplanted pinecones is not relevant or where we think 
the context makes it clear which sort of pinecone we intend, we omit the modifier and simply 
use the word "pinecone". 
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Figure 4. Some pinecones of width 15. a. The dashed edges may belong to 
the graph, or not. b. A pinecone. c. A closed pinecone. 

Figures [Ho and [He show two specific ways to make the choices indicated in Figure [Ha and 
obtain a pinecone of width 15. The pinecone of Figure [He is closed, meaning that it contains 
no vertex of degree 1. (Such a vertex can only occur at the right border of the pinecone, and 
occurs if the rightmost vertex of some horizontal segment does not belong to a vertical edge, as 
shown in Figure [Hp-) An Aztec diamond graph is an example of a closed pinecone. Let us color 
the cells of the square lattice black and white in such a way that the cell containing the vertices 
(0, 0) and (1, 1) is black. The faces of the pinecone are the finite connected components of the 
complement of the graph in R 2 . The faces of a pinecone P are of three types: black squares, 
white squares, and horizontal rectangles consisting of a black cell to the left and a white cell 
to the right. We insist on the distinction between a cell (of the underlying square lattice) and 
a square (a face of P that has 4 edges). For instance, the longest row of a pinecone of width 
2k — 1 contains exactly 2k — 1 cells, but may contain no square at all. Denoting by i (resp. r) 
the leftmost (resp. rightmost) cell of the longest row of P, we say that P is rooted on (£, r). (If 
P is standard, then I is the cell with (0,0) as its lower-left corner.) 

It is easy to see that a pinecone P is closed if and only if the rightmost finite face of each 
row is a black square. In this case, the rightmost black square in each row is also the rightmost 
cell of the row. If moreover P is standard, it is completely determined by the position of its 
black squares. Equivalently, it is completely determined by the position of its odd vertical edges. 
Conversely, consider any finite set S of black squares whose lower-left vertices lie in the 90 degree 
wedge bounded by the rays y = x > and — y = x > 0. Assume that S is monotone, in the 
following sense: the rows that contain at least one square of S are consecutive (say from row —j 
to row i) and for m > (resp. m < 0), the rightmost black square in row m occurs to the left 
of the rightmost black square in row m — 1 (resp. m + 1). Then is a (unique) closed standard 
pinecone whose set of black squares is S. 

We shall often consider the empty graph as a particular closed pinecone (associated with the 
empty set of black squares). The empty graph has one perfect matching, of weight 1. 

2.4. The core of a pinecone 

When a pinecone P is not closed, some of the edges of P cannot belong to any perfect matching 
of P. Specifically, if v is a vertex of degree 1 in P, then in any perfect matching of P, v must 
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be matched with the vertex to its left (call it u) , so that u cannot be matched with any of its 
other neighbors. Indeed, there can be a chain reaction whereby a forced edge, in causing other 
edges to be forbidden, leads to new vertices of degree 1, continuing the process of forcing and 
forbidding other edges. An example of this is shown in Figure [5l The left half of the picture 
shows a non-closed pinecone P, and the right half of the picture shows a closed sub-pinecone 
P of P along with a set of isolated edges. The reader can check (starting from the rightmost 
frontier of P and working systematically leftward) that each of the isolated edges is a forced 
edge (that is, it must be contained in every perfect matching of P), so that a perfect matching 
of P is nothing other than a perfect matching of P together with the set of isolated edges shown 
at right. In this subsection, we will give a systematic way of reducing a pinecone P to a smaller 
closed pinecone by pruning away some forced and forbidden edges. 



EZI 



Figure 5. From a pinecone P to its core P. 

It can easily be checked that the union or intersection of two standard pinecones is a standard 
pinecone, and that the union or intersection of two closed standard pinecones is a closed standard 
pinecone. It follows that, if P is a standard pinecone, there exists a largest closed standard sub- 
pinecone of P, namely, the union of all the closed standard sub-pinecones of P. We call this 
the core of P and denote it by P. (If P is not a standard pinecone but a transplanted pinecone 
rooted at the cell with lower-left corner (a, b), we define Pq as the standard pinecone obtained 
by translating P by (—a, —b), and we define the core of P as the core of Pq translated by (a, b). 
However, for the rest of this section we will restrict attention to standard pinecones.) 

Here is an alternative (more constructive and less abstract) approach to defining the core. 
Let P be a standard pinecone. Let b be the rightmost black square in row of P, let b\ be the 
rightmost black square in row 1 of P that lies strictly to the left of bo, let &2 be the rightmost 
black square in row 2 of P that lies strictly to the left of b±, and so on (proceeding upwards); 
likewise, let 6_i be the rightmost black square in row — 1 of P that lies strictly to the left of bo, 
and so on (proceeding downwards). If at some point there is no black square that satisfies the 
requirement, we leave b m undefined. Consider all the faces of P that lie in the same row as, and 
lie weakly to the left of, one of one of the bk's. This set of faces gives a closed pinecone P. At 
the same time, it is clear that any closed sub-pinecone Q of P must be a sub-pinecone of P. For, 
the rightmost black square in row of Q can be no farther to the right than bo, which implies 
that the rightmost black square in row 1 of Q can be no farther to the right than b\, etc.; and 
likewise for the bottom half of Q. Hence the sub-pinecone P we have constructed is none other 
than the core of P as defined above. 

If P is closed, then P = P. Note that a closed pinecone always admits two particularly simple 
perfect matchings: one consisting entirely of horizontal edges, and the other consisting of the 
leftmost and rightmost vertical edges in each row (and no other vertical edges) along with some 
horizontal edges (Figure EJ)- In particular, the rightmost vertical edges of a closed pinecone are 
never forced nor forbidden. 

Let P be a pinecone with core P. There is a unique perfect matching of P \ P consisting 
exclusively of horizontal edges (see Figure [5]) ; let H be the edge set of this perfect matching. 
Every perfect matching of P can be extended to a perfect matching of P by adjoining the edges 
in H, so m(P) > m(P). We now show that every perfect matching of P is obtained from a 
perfect matching of P in this way. 
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Figure 6. Two particularly simple matchings of a closed pinecone. 
Proposition 3. Let P be a pinecone with core P. Then m{P) = m(P). 

Proof. We will prove this claim by using a procedure that reduces a sub-pinecone Q of P to a 
smaller sub-pinecone with the same matching number. Let Q be a sub-pinecone of P whose core 
coincides with P. If Q is not closed, then there must be at least one vertex of degree 1 along the 
right boundary of Q. Let v = (a, b) be one of the the rightmost vertices of degree 1 in Q. Then 
v is the rightmost vertex in one of the rows of Q. Assume for the moment that v lies strictly 
above the longest row of Q (that is, b > 1). See the top part of Figure[7]for an illustration of the 
following argument. Let u be the vertex to the left of v. Then the edge joining u and v is forced 
to belong to every perfect matching of Q, while every other edge containing u is forbidden from 
belonging to any perfect matching of Q. Hence the graph Q' obtained from Q by deleting u, v, 
and every edge incident with u or v has the same matching number as Q. Furthermore, Q' is a 
pinecone, unless the vertex V\ = (a — 1, b + 1) belongs to Q. In this case, v\ has degree 1 in Q' . 
Let i be the largest integer such that Vj = (a — j,b + j) belongs to Q for all < j < i. Applying 
the deletion procedure to the vertices v — vq, v±, . . . , (in this order) yields a pinecone Q*. 
Assume now that 6=1. Applying the deletion procedure to all the vertices of Q of the form 
(a — j, 1 + j) or (a — j, — j) yields again a pinecone Q* (see Figure bottom). By symmetry, we 
have covered all possible values of b. 




Figure 7. Some sequences of edge-deletions starting and ending with a pinecone. 

Observe that m(Q) = m(Q*). Additionally, we can check that the core of Q* is P. The 
only thing we might worry about is that in passing from Q to Q* , we removed some edges that 
belong to P. The examination of Figure [5] shows that we would have, in particular, removed the 
rightmost vertical edge is some row of P. However, this cannot happen, because the removed 
edges were all forced or forbidden, whereas the rightmost edges of P are neither forced nor 
forbidden (Figure [6|). 

To prove that m(P) = m(P), take Q = P and use the preceding operation repeatedly to con- 
struct successively smaller graphs Q* , Q** , . . .such that m{P) = m(Q) = m(Q*) = m(Q**) = 
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• • ■ and P = Q = Q* = Q** = ■ ■ ■ . Eventually we arrive at a closed sub-pinecone of P whose 
core is P; that is, we arrive at P itself. And since each step of our construction preserves m(Q), 
we conclude that m(P) = m(P), as claimed. □ 

2.5. A CONDENSATION THEOREM FOR CLOSED PINECONES 

Let P be a closed pinecone, with longest row consisting of 2n + 1 squares. Let A be the 
smallest diamond graph containing P (the longest row of A contains exactly 2n + 1 cells) . Let 
G denote the spanning subgraph of A whose edge-set consists of all edges of P, all horizontal 
edges of A, and all even vertical edges of A (Figure [8]). Observe that G is a pinecone. Moreover, 
among the spanning subgraphs of A that are pinecones and contain P, G has strictly fewer edges 
than the others. Since no odd vertical edge is added, P is actually the core of the pinecone G. 



P 



G 



I 



I 



Figure 8. Completing a pinecone P into a spanning pinecone of an Aztec 
diamond graph. 

Let us now use the notation of Corollary [2 That is, Gn = G fl An, and so on. Then 
Gn,Gs,Gw,Ge and Gc are (standard or transplanted) pinecones. Let P , P s , P w , P E and 
P c denote their respective cores. (These are not to be confused with Pjv, etc., which are the 
intersections of P with A N , etc.) We will often call P N , P s , P w , P E and P c "the five sub- 
pinecones" of P, even though, strictly speaking, P admits other sub-pinecones. An example is 
given in Figure [9l Let £q (resp. ro) be the leftmost (resp. rightmost) cell of the longest row i?o 
of P. Similarly, let n (resp. r_i) denote the rightmost cell of the row just above (resp. below) 
Ro- Observe that the cells r ,ri and r_i correspond to black squares of P. Finally, let £' Q be 
the black cell of Rq following £ , and let r' the black square of i?o preceding ro (if it exists). 
In light of the basic properties of the core (both the abstract definition and the algorithmic 
construction), we can give the following alternative description of the five sub-pinecones of P. 

Proposition 4. Let P be a closed pinecone. With the above notation, P N (resp. P ) is the 
largest closed sub-pinecone of P whose rightmost cell is n (resp. r-i). Similarly, P w (resp. P E ) 
is the largest closed sub-pinecone whose rightmost (resp. leftmost) cell is ?* (resp. t' ). Finally, 
P c is the largest closed sub-pinecone rooted on (£'o,r' Q ). 

This proposition implies that a pinecone P that is neither empty, nor reduced to a black 
square can be reconstructed from its four main sub-pinecones P N , P s , P E and P w . Indeed, 
the part of P located strictly above its longest row coincides with the top part of P . More 
precisely, row r of P, with r > coincides with row r — 1 of P N . Similarly, for r < 0, row r 
of P coincides with row r + 1 of P s . It thus remains to determine the longest row of P. This 
row is obtained by adding a 2-by-l rectangl^l to the left of the longest row of P E , and then 
superimposing the longest row of P w . 



2 This rectangle is actually only useful if P w is empty or reduced to a single black square. 
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Figure 9. The five sub-pinecones of a pinecone P. 



Let us now apply Corollary [2] to the graph G obtained by completing P into a spanning 
pinecone of A (Figure IBJl- By Proposition [31 since P is the core of G, m{G) = m(P), and similar 
identities relate the matching numbers of Gn and P N , etc. 

Theorem 5 (Condensation for closed pinecones). The matching number of a closed 
pinecone P is related to the matching number of its closed sub-pinecones by 

m(P)m{P c ) = m{P w )m(P E ) + m{P N )m{P s ). 

We will state in Section [4] a more general condensation result dealing with the matching poly- 
nomial, rather than the matchings number, of closed pinecones (Theorem [T3jl . 

3. Pinecones for the Gale-Robinson sequences 

The pinecones introduced in the previous section generalize Aztec diamond graphs. The 
number of perfect matchings of the diamond graph of width 2n — 3 is the nth term in the 
recurrence 

a(n)a(n — 2) = a(n — l)a(n — 1) + a(n — l)a(n — 1), 

with initial conditions a(0) = a(l) = 1. More generally, the three-term Gale- Robinson sequences 
are governed by recurrences of the form 

a(n)a(n — m) = a(n — i)a(n — j) + a(n — k)a(n — £), (1) 

with initial conditions a(n) = 1 for n = 0, 1, . . . ,m — 1. Here, i,j, k and £ are positive integers 
such that i + j = k + £ = m, and we adopt the following (important) convention 

j = mm{i,j,k,£}. 

Our purpose in this section is to construct a sequence of (closed) pinecones (P(n)) n >o = 
(P(n;i,j,k,£)) n>0 for each set of parameters {i,j,k,£} such that i + j = k + £ = m, and to 
show that the matching numbers of the pinecones in our sequence satisfy the corresponding 
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Gale-Robinson recurrence. More specifically, our family of graphs will be constructed in such a 
way that 

• P(n) c is P(n — m) transplanted to (2, 0) (that is, shifted two steps to the right), 

• P(n) w is P{n-i), 

• P(n) E is P(n — j) transplanted to (2,0), 

• P(n) N is P(n — k) transplanted to (1, 1), and 

• P(n) s is P(n — £) transplanted to (1, —1). 

In our construction, we use the fact that a closed pinecone is completely determined by its set 
of odd vertical edges, that is, vertical edges of the form {(a, 6), (a, b + 1)} where a + b is odd. 
We introduce two functions, an upper function U and a lower function L, which will be used 
to determine the positions of the odd vertical edges in the (closed) pinecone P(n;i,j,k,£): for 
r > 0, let 

mc + kr + i — n — 1 

U(n,r,c) = 2c + r-3-2 

J 

(2) 

mc + £r + i — n — 1 

L(n,r,c) = 2c + r-3-2 

3 

Observe that the parameters k and £ play symmetric roles. Also, U(n, 0, c) = L(n,0, c). The 
function U will describe the upper part of the pinecone, while L will describe its lower part. 
Recall that, by convention, the longest row of a standard pinecone is row and its South- West 
corner lies at coordinates (0,0), as shown in Figure 21 

To locate the vertical odd edges in row r > 0, calculate the values U(n, r, c) for c = 0,1,... 
This will be a (strictly) decreasing sequence, since m > 2j (recall that i + j = m and j < i). 
Retain those values U(n,r,c) that are larger than r, and place a vertical edge in the rth row 
at abscissa U(n,r,c), that is, an edge connecting (U(n,r,c),r) and (U(n, r, c),r + 1) (an odd 
edge, since U(n, r,c) + r is odd). The first row not containing such an edge (and therefore not 
included in the pinecone) is the first one for which U(n, r, 0) < r. Observe that U(n, r, 0) — r 
is a decreasing function of r (since j < k). This property guarantees that if the rth row is 
empty, then all higher rows are empty too. It also implies that the rightmost vertical edge in 
row r (which is located at abscissa U(n, r, 0)) lies to the right of the rightmost vertical edge in 
row r + 1. (To see this, note that U(n, r, 0) — r is always an odd number. So the inequality 
£7(n, r + 1, 0) - (r + 1) < U(n, r, 0) - r implies U(n, r + 1, 0) - (r + 1) < U(n, r, 0) - r - 2, or 
U(n,r + 1,0) < U(n, r, 0). That is, the set of odd edges (equivalently, of black squares) given 
by Formula ^ satisfies the "top part" of the monotonicity condition described at the end of 
Section 12.31 the rightmost odd edge in row r > 0, if it exists, lies to the left of the rightmost 
odd edge in row r — 1. 

Similarly, to locate the edges in row — r < 0, calculate the values of L(n, r, 0) > L(n, r, 1) > • • ■ 
and retain those larger than r. For each, place a vertical edge in row — r at abscissa L(n,r,c), 
that is, connecting (L(n, r, c), — r) and (L(n,r,c),—r + 1). Observe that L(n, 0, c) = U(n, 0,c), 
so that the collection of odd vertical edges in row is the same whether it is determined from 
U or from L. 

The monotonicity properties satisfied by the positions of the odd edges imply that there exists 
a unique standard closed pinecone whose set of odd vertical edges coincides with the set we have 
constructed via the functions U and L. To obtain this pinecone, draw horizontal edges from 
(r, r + 1) to (U(n, r, 0), r+ 1) and from (— r, — r) to (L(n, r, 0), — r) for all r > 0. Finally, place all 
the appropriate even vertical edges. Since these steps are so routine, we regard the pinecone as 
fully described once the set of odd vertical edges has been specified. This point of view simplifies 
the exposition. 

Observe that P(n) is empty if and only if U(n, 0, 0) < 0, which is equivalent to U(n, 0, 0) < —1 
(since U(n, 0, 0) is odd), which is easily seen to be equivalent to n < m (using the fact that 
m = i + j). 



PERFECT MATCHINGS FOR THE THREE-TERM GALE-ROBINSON SEQUENCES 



13 



Example. Take (i,j,k,£) = (5,2,3,4) and determine P(12). The above definition of U and L 
specializes to 



U(n, r, c) 
L(n,r, c) 



2c + r - 3 - 2 
2c + r - 3 - 2 



7c + 3r - 8 



7c + 4r - 8 



In row 0, we find odd edges with lower vertices (5, 0) and (1, 0). In row 1, there is one odd edge 
at (4,1). This is the top row of the diagram because £7(12,2,0) = 1 < 2. Turning to the lower 
portion of the diagram, there is one odd edge with lower vertex (2, —1) and none in row —2 or 
below. Completing the diagram is now routine, and gives the pinecone P(12) which is shown 
in Figure flOl together with its 14 perfect matchings. Accordingly, the Gale- Robinson sequence 
a(n) associated with (5,2,3,4) satisfies a(12) = 14. 
A larger example is presented after Corollary [TOl 



<P« 

(2) (2) (4) 



(0,0)' 



Figure 10. The pinecone P(12; 5, 2, 3, 4), with black squares indicated, and its 
14 perfect matchings (a cross stands for any of the two matchings of a square) . 




The pinecones based on the functions U and L satisfy a remarkable property: the odd edges 
(or, equivalently, the black squares) in rows r and r + 1 are interleaved. That is, between two 
black squares in row r > 0, there is a black square in row r + 1, and similarly, between two black 
squares in row r < 0, there is a black square in row r — 1. This can be checked on the small 
example of Figure [TUl but is more visible on the bigger example of Figure fT2l 

Lemma 6 (The interleaving property). For all values of n,r and c, the functions U and L 
defined by (j2]) satisfy 



and 



Proof. We have 



U(n, r, c+ 1) + 1 < U(n, r + 1, c) < U(n, r, c) - 1 



L(n, r, c+ 1) + 1 < L(n, r + 1, c) < L(n, r, c) — 1. 



U(n,r+l,c)-U(n,r,c+l) = 2 
But 



mc -\-kr + i — n — 1 + m 



-2 



mc + kr + i — n — 1 + k 



-1. 



mc + kr + i — n — 1 + m mc + kr + i — n — 1 + k 



> 1, 



3 3 3 

so that the two floors occurring in the above identity differ by 1 at least. Consequently, 

U(n, r + 1, c) - U{n, r, c+ 1) > 2 - 1 = 1. 

The three other inequalities are proved in a similar manner. 



□ 
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We now wish to apply the condensation theorem (Theorem [5]) to the pinecones P(n) we have 
just defined. Using the notation of Theorem [H we will verify that, up to translation, P(n) w = 
P(n-i), P(n) E = P(n-j), P(n) N = P(n-fc), P(n) s = P(n-£) and P(n) c = P(n-m). These 
equivalences will follow from the interleaving property and the following algebraic equalities. 

Lemma 7. For any choice of parameters (i,j,k,£), the functions U and L defined by (J2j) satisfy: 

U(n — i,r,c— 1) = U(n,r,c), L(n — i,r,c—l) = L(n,r,c), 

U{n — j,r,c) = U(n,r,c)—2, L(n — j,r,c) = L(n,r,c) — 2, 

U(n — k,r—l,c) = U(n,r,c) — 1, L(n — £, r — 1, c) = L(n, r, c) — 1, 

U(n — £, r + 1, c — 1) = U(n, r, c) — 1, L(n— k,r + l,c— 1) = r, c) — 1, 

U(n — m,r,c—l) — U(n,r,c)—2, L(n — m,r,c—l) — L(n,r,c) — 2. 

Proof. The L-identities are symmetric to the [/-identities upon exchanging k and i, so that there 
are really 5 identities to prove. These can all be verified by routine algebraic manipulations. Let 
us check for instance the fourth identity satisfied by U: 

' m(c - 1) + k(r + 1) + % - (n -£)-l 



U(n-£,r + l,c-l) 



2(c-l) + (V + l)-3-2 

mc + kr + i 

2c + r 



J 

l-m + k + £ 



2c + r-4-2 
U(n, r, c) — 1. 



J 

mc + kr + i — n — 1 
3 



since to 



k+£ 



We leave it to the reader to verify the remaining 4 identities. 



□ 



We now check that these identities imply that the pinecones are related to one another as 
claimed. 

Proposition 8. Let P(n) = P(n:i,j,k,£) be the sequence of pinecones associated with the 
parameters (i, j, k, £). Then for n > m, the five closed sub-pinecones of P(n) satisfy 

P(n) w = P(n-i), P{n) E = P{n~ j), 

P(n) N = P(n- k), P(n) s = P(n - £), 



P{n 



and 

P(nf 

These identities hold up to a translation. 

Proof. Begin by checking that P(n) w = P(n — i). Using the description of P(n) w given in 
Proposition IH and the fact that the black squares of P(n) are interleaved, we see that the odd 
vertical edges in P(n) w are those of P{n), except that the rightmost odd edge in each row 
has been removed. (If this was the only odd edge in the row, then the entire row disappears.) 
Therefore P(n) w can be constructed by following the construction for P(n), but beginning with 

c = 1,2,... instead of c = 0, 1, This means that in row r > of P(n) w , odd edges appear 

at positions U(n, r, 1), U(n, r, 2), . . . as long as these values continue to exceed r. (Similarly in 
rows r < 0, using L instead of U.) 

Let us now compare this with P(n — i). In row r > of P(n — i) } odd edges appear in positions 
U (n — i, r, 0), U(n — i, r, 1), . . ., as long as these values continue to exceed r. However we showed 
in Lemma[7]that U(n — i,r,c — 1) = U(n,r,c), so the sequence of odd edges in row r is the same 
in P(n) w and in P(n — i). The situation is similar in rows r < using the equality for L. As 
we remarked above, a pinecone is determined by its odd edges (and the position of its leftmost 
edge), so P(n) w = P(n — i). 

The other four equivalences are similar. The only new development is that instead of being 
positioned at the origin, the smaller pinecones are now offset by one or two columns (in all 
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four cases) and possibly rows (in the case of P{n) N and P(n) s ). We will look at P(n) s as an 
example, and let the reader supply the details for the remaining three cases. 

As noted after Proposition [4j for r < 0, row r of P(n) s coincides with row r — 1 of Pin) 
(see Figure [9] for an example). For r > 0, the leftmost cell of row r of P{n) s lies two steps to 
the right of the leftmost cell of row r — 1 of P(n). Moreover, the interleaving property implies, 
by induction on r, that the last (i.e., rightmost) black square of row r of P(n) c is the next-to 
last black square of row r — 1 of P(n). Thus the odd edges of P(n) s are located as follows: for 
rows — ?*, with r = 1, 2, . . ., in columns L(n, r, 0), L(n, r, 1), . . ., as long as these numbers continue 
to exceed r; and for rows r — 0, 1, 2, . . ., in columns U (n, r, 1), U(n, r, 2), . . ., as long as these 
numbers continue to exceed r + 2. 

Let us now look at a copy of P(n — t) positioned with its origin at (1,-1). After this 
translation, the odd vertical edges in rows — r, with r = 1,2,... are located at abscissas L(n — 
£,r — l,c) + 1, for c > and as long as these numbers continue to exceed r. Lemma [7| then 
implies that the bottom parts of P(n) s and of the translate of P(n — £) coincide. After the 
translation, the odd vertical edges of P(n — t) lying in rows r, with r = 0, 1, 2, . . . are located at 
abscissas U(n — I, r + 1, c) + 1, for c > and as long as these numbers continue to exceed r + 2. 
Lemma [7] then implies that the top parts of P(n) s and of the translate of P(n — I) coincide. 

This completes the analysis for P(n) s ; the verifications for P(n) N , P(n) E and P(n) c are 
similar (and even identical, up to symmetry, in the case of P(n) N ). □ 



P(4) = □ 



P(5) = I ■ 
P(5) E = P(4) 



P(6) = C 

P(6) s = P(5) 



P(7) = 



P(8) = 



P(9) = C 



I 



P(7) w =P(4), P(7) s = P(6) 
I 



P(10) = 



T 



P(8) JV = P(8) s = P(6) 



P(9) JV = P(7) 



I 



I 



P(10) w ' = P(7) 



Figure 11. Recursive graphical construction of the pinecones associated with 
the Somos-4 sequence. At each stage, one (or two) of the components that are 
superimposed to form the pinecone is highlighted. 

Remark: a recursive construction of the pinecones P(n). The above proposition, com- 
bined with Proposition HI provides an alternative way of constructing the sequence of pinecones 
P(n) associated with a given set of parameters (i,j,k,£). For < n < m, we put P(n) equal 
to the empty graph (which has one perfect matching), and for m < n < m + j, we put P(n) 
equal to the graph with four vertices and four edges surrounding one square face (which has 2 
perfect matchings). Then, for n > m + j, it suffices to superimpose P(n — i), P(n — j), P(n — k) 
and P(n — £), and add a 2-by-l rectangle to the left of the longest row of P(n — j). More 
precisely, the four above pinecones must be positioned in such a way the leftmost cell of P(n — i) 
(resp. P{n-j), P(n-k), P(n-l)) has its South- West corner at (0,0) (resp. (2,0), (1,1), (1,-1)), 
while the 2-by-l rectangle has its South- West corner at (0,0). (This rectangle is only needed if 
P(n — i) is empty or consists of a single black square. Typically this 2-by-l rectangle comes for 
free as part of P[n — i). Note that we do not claim that this rectangle is a face of the pinecone; 
the odd edge joining (1, 0) and (1, 1) will be present or absent in P(n), according to whether it 
is present or absent in P(n — i).) This gives a graphical, inductive way of constructing P(n). 
This method is illustrated in Figure [TT] by the case of the Somos-4 sequence, for which 

a(n)a(n — 4) = a(n — 3)a(n — 1) + a(n — 2) 2 . 
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That is, (i, j, k, £) = (3, 1, 2, 2) and m = 4. 

We can now state our combinatorial interpretation of the Gale-Robinson numbers. 

Theorem 9. Let P(n) = P{n;i, j,k,i) be the sequence of pinecones associated with the param- 
eters (i,j,k,£). Let a{n) denote the number of perfect matchings of P(n). Then a(n) = 1 for 
n < m, and for n > m, the sequence a(n) satisfies the following Gale-Robinson recurrence: 

a(n)a(n — to) = a(n — i)a(n — j) + a(n — k)a(n — £). 

Proof. We have already observed that the pinecone P(n) is empty for n < m. Hence the initial 
conditions apply correctly. Now for n > m, Theorem \E\ states that the matching matching of 
P(n) is related to the matching numbers of its closed sub-pinecones by 

m(P{n))m(P{nf) = m(P{n) w )m(P{n) E ) + m(P{n) N )m{P{n) s ). 



Proposition [8] then implies that m(P(? 



\c\ - 



m(P(n — to)), etc. Therefore, 



□ 



m(P(n))m(P(n — to)) ~ m(P(n — i))m(P(n — j)) + xm(P(n — k))m(P(n — £)), 

which is the recurrence relation satisfied by a{n). 

Before we study a specific example, let us state an obvious corollary of Theorem [9l 

Corollary 10. Let i,j,k,£ be positive integers such that i + j = k + £ = m. The recurrence 
relation 

a(n)a(n — to) = a(n — i)a(n — j) + a(n — k)a(n — t), 
with initial conditions a(n) = 1 for n < m, defines a sequence of positive integers 

Example. We give a specific example in the case where (i,j,k,£) = (6,2,5,3) and 
n = 25. We also show how to use the VAXmaple software package (available at 



http://jamespropp.org/vaxmaple.cl to compute the number of perfect matchings in the con- 



structed pinecone, which can be seen to be the 25th term in the appropriate Gale-Robinson 
sequence. 

Considering first the upper portion of P{n), we fix r and then consider the first few values of 
U(n, r, c) as c = 0, 1, 2,3,.. .: 

r = 0: {17,11,5,-1,...} 
r=l: {14,8,2,-4,...} 
r = 2: {9,3,-3,-9,...} 
r = 3: {6,0,-6-12,...} 
r = 4: {1,-5,-11,-17,...} 

Since the c = value for r = 4 is already less than r, there are only three non-empty rows above 

the middle (longest) row in this pinecone. For the lower portion of the diagram, we obtain the 
following values of L(n, r, c): 

0: {17,11,5,-1,...} 

1: {16,10,4,-4,...} 

2: {13,7,1,-5,...} 

3: {12,6,0,-6,...} 

4: {9,3,-3,-9,...} 

5: {8,2,-4,-10,...} 

6: {5,-1,-7,-13,...} 

Completing the construction, we arrive at the graph of Figure [121 

It is easy to translate this into the format required by the computer program VAXmaple, 
written by Greg Kuperberg, Jim Propp and David Wilson to count perfect matchings of finite 
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L_U 

Figure 12. The pinecone P(25; 6, 2, 5, 3). 

subgraphs of the infinite square grid. In this format, each vertex present in the graph is repre- 
sented by a letter. The choice of letter indicates whether any edges are omitted when connecting 
the vertex to its nearest neighbours — each vertex having up to four of these. An X indicates that 
no edges are omitted; an A indicates that the edge leading upward from the vertex is omitted; 
a V indicates the omission of the downward edge. (For a more detailed explanation of the soft- 
ware, see http://jamespropp.org/vaxmaple.doc.) The encoding of the pinecone of Figure fT2l 
is given in Figure fl3l 

XVXX 
XXAVXVXX 
XXXVAVAXXVXVXX 
XVXVAXAVXVAXAVXVXX 
XAVAXXVAVAXXVAVAXX 
XAVXVAXAVXVAXAXX 
XAVAXXVAVAXX 
XAVXVAXVXX 
XAVAXX 
XAXX 

Figure 13. The pinecone P(25; 6, 2, 5, 3) as a VAX file. 

Counting the perfect matchings in this pinecone by running the above input through the 
VAXmaple program and then through Maple produces 167,741, as it should, since the 25th term 
of the Gale-Robinson sequence constructed from (6,2,5,3) is 167,741. 

4. The Gale-Robinson bivariate polynomials 

As stated in Corollary [lOl Theorem [9] implies that the three-term Gale- Robinson sequences 
consist of integers. In this section, we refine this result as follows. 

Theorem 11. Let i, j, k, £ and m be positive integers such that i + j = k + £ = m. Let w and z 
be two indeterminates, and define a sequence p(n) = p(n; w, z) by p{n) = 1 for n < m and for 

n > m, 

p{n)p{n — m) = wp(n — i)p(n — j) + zp(n — k)p(n — £). 
Then p(n) is a polynomial in w and z with nonnegative integer coefficients. 

The proof goes as follows: we have already seen that p(n; 1, 1) counts perfect matchings 
of the pinecone P = P(n;i,j,k,£) constructed in Section 03 We will prove that p(n;u 2 ,v 2 ) 
counts these matchings according to two parameters. More precisely, we begin by giving in 
Section 14.11 a condensation theorem that computes inductively the matching polynomial (rather 
than the matching number) of closed pinecones. We observe that this theorem takes a simpler 
form when applied to interleaved pinecones (a class of pinecones that contains all Gale- Robinson 
pinecones). In Section [OJ we define the special horizontal edges of a pinecone. We then define 
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the partial matching polynomial of a pinecone P as the matching polynomial M(P) in which the 
weights of non-special edges are set to 1. We specialize the condensation theorem of Section |4~T1 
to a condensation theorem for the partial matching polynomial of interleaved pinecones. Its 
application to the Gale- Robinson pinecones P(n;i,j,k,£) implies that the polynomial q(n) = 
q(n; u, v) that counts perfect matchings of P(n) according to the number of vertical edges (the 
exponent of v) and horizontal special edges (the exponent of u) satisfies q(n) = 1 for n < m and 

q{n)q(n — m) = u qin — i)q(n — j) + v q(n — k)q(n — i) 

for n > m. This shows that q(n; u, v) = p(n; u 2 , v 2 ) and implies Theorem fTTl 

4.1. A CONDENSATION THEOREM FOR THE MATCHING POLYNOMIAL 

Let us go back to the condensation theorem for closed pinecones (Theorem [5]) . We now state 
and prove a stronger result dealing with the matching polynomial rather than the matching 
number. Let P be a closed pinecone and A the smallest diamond graph that contains it, with 
G defined as in the beginning of Section l2~5l and with n, s, e, w as in Theorem Q] and Figure [H 
Since P is the core of G, the matching polynomial M{G) equals M(P)M(G\P). Similar results 
hold for the sub-pinecones P c , P w , P E , P N and P s . Corollary [2] gives: 

M(P)M(G \ P)M(P C )M(G C \ P C ) = nsM(P w )M(G w \ P W )M(P E )M{G E \ P E ) 

+ ewM(P N )M(G N \ P N )M(P S )M(G S \ P s ). (3) 

Since P is the core of G, the graph G \ P has a unique perfect matching, which is formed of 
horizontal edges only. Hence M(G \ P) is a monomial. The same holds for the other graph 
differences occurring in (f3]). We can thus rewrite this identity as 

M(P)M(P C ) = aM(P w )M{P E ) + /3M(P N )M(P S ) 

for some Laurent monomials a and j3 (indeed, negative exponents may arise from the division by 
M(G \ P)M(Gc \ P C ))- Our objective in this subsection is to prove that these monomials only 
involve nonnegative exponents (so that they are ordinary monomials), and to describe them in 
a more concise way. 

We introduce the following definition, illustrated in Figure HH 

Definition 12. Let P be a closed pinecone. A horizontal edge is a left edge if it is the leftmost 
horizontal edge in the horizontal segment of P that contains it. 

A horizontal edge with leftmost vertex is even ifi + j is even, odd otherwise. 



Figure 14. The left edges of a closed pinecone. 



For a (standard) pinecone P, we denote by P> (resp. P>) the pinecone formed of rows 
0, 1, 2, . . . of P (resp. rows 1,2,.. .). We use similar notations for the bottom part of P. These 
definitions are extended to transplanted pinecones in a natural way: for instance, (P N )>, which 
we simply denote P> , consists of rows 1, 2, ... of P. 

Observe that, for any pinecone P, 

P^cPf while P^CP^. (4) 
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Figure 15. The inclusion properties P< C P< and P> C P> (the small 
subpinecones are dashed). The leftmost figure also shows some edges of the 
horizontal matching of P< \ P< . 



Both properties are illustrated in Figure H21 Consequently, the graph difference P c \ P N is 
formed of edges that lie in P< , and P c \ P N C P< \ P< . Let us describe more precisely the 
horizontal edges of these two graph differences. For j < 1, if there are any horizontal edges of 
P< \ P< ly m S a t ordinate j, then the number of them is odd, say 2k j + 1, and these edges are 
the 2k j + 1 rightmost horizontal edges of P c found at ordinate j (Figure [15j left). If j < 0, 
all these edges belong to P c \ P N . However, for j = 1, only a subset of these edges, of even 
cardinality, belong to P C '\P N . (In the example of FigurefTBl the two leftmost thick edges shown 
at ordinate 1 do not belong to P c \ P N .) The graph P< \ P< thus has a unique horizontal 
matching, which has kj + 1 edges at ordinate j < 1. We denote by H~ (P c \ P N ) the product of 
the edges of this matching having ordinate < 0. The fact that the horizontal edges of P c \ P N 
found at ordinate j coincide with those of P< \ P^f allows us to use the notation H~(P C \ P N ) 
rather than something like H~(P< \ P^) which would have been heavier. 
Symmetrically, 

P| C P§ while P? C P|, (5) 

so that the graph P C \P S lies in P> . We denote by H + (P C \P S ) the product of the edge- weights 
of the horizontal matching of P> \ P> lying at a positive ordinate. 



rightmost cell of P 




rightmost cell of P s 

Figure 16. Left: The edges occurring in the first term of the refined conden- 
sation theorem, with the pinecone P w shown. Right: The edges occurring in 
the second term. Here H + (P C \ P s ) = 1. The two distinguished pinecones are 
P c and P N (dashed). 

We can now state a condensation theorem for the matching polynomial of closed pinecones. 
See Figure [16] for an illustration. 
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Theorem 13 (The matching polynomial of closed pinecones). The matching polynomial 
of a closed pinecone P is related to the matching polynomial of its sub-pinecones by 



M(P)M (P c ) = I ' [ a)M (P w ) M (P E ) 



( II a ) H-(P C \P N )H+(P C \P S )M (P N ) M (P s ) 



where 

• E\ is the set of left edges of P not belonging to P w , 



Ei is the union of three edge-sets E\, for < i < 2; 

— i?" = {e,w} contains the eastmost and westmost vertical edges of P, 

— E\ contains the even edges at ordinate 1 not belonging to P N , 

— E\ contains the odd edges at ordinate not belonging to P s . 



Proof. In this proof, we adopt the following notation: for each edge set E, we also denote by E 
the product of the edges of the set. For a graph G having a unique horizontal (perfect) matching, 
we denote this matching by H{G). 

Let us return to (J3j) . Recall that G\P has a unique matching, consisting of horizontal edges 
only. Denoting by A the smallest diamond graph containing P, we observe that M(G \ P) = 
H(A)/H(P) (see Figure [8|. Similar identities hold for the other pinecones occurring in ([3]). For 
instance, M(G W \ P w ) = H(A W )/H(P W ). This allows us to rewrite © as 

M(P)M(P°) = ns H ^ H ^ H{P)H(PC) M{pW)M{pE) 
[ ' [ ' H{A)H{A C ) H{P W )H(P E ) [ ' { ' 

H{A„)H(A S ) H(P)H(PC) M{pN)M{p s y 
H(A)H(A C ) H(P N )H(P S ) y ' v ' 

Let us begin with the two factors involving A and its subgraphs. It is easy to see, with the help 
of Figure H, that 

H(A W )H(A E ) _ 
US H(A)H(A C ) L 

The second factor involving A, namely ewH (An)H (As) / (H (A)H (Ac)), is a multiple of e and 
w (all the other edges are horizontal) and thus cannot be equal to 1. Denoting by the graph 
formed by the even horizontal edges lying at ordinate 1, and introducing similar notations 



to r.W v 

(e) r (o) 



Lq and Lq' , one finds 



H(A N )H(A S ) L\ L { 
H(A)H(A C ) L (o) L (e)- ^ 

It remains to describe the two factors that involve P and its subgraphs. For the first one, we note 
that H(P)/H(P E ) is simply the product of the left edges of P. Similarly, as P c = (P W ) E , the 
ratio H(P W )/H(P C ) is the product of the left edges of P w . This gives the following expression 
for the first factor: 

H(P)H(P C ) n 
H(P W )H(P E ) 11 ' 

with Ei defined as in the theorem. 

To express the second factor involving P, let us first separate in H(P) the edges that lie at 
ordinate j = 0, j = 1, j > 1, j < 0. This gives 



H(P) = L e) • L[ o) 



H(P?) H(Pi) 



L[ e) nP N L { o) nP s ' 
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For the other 3 pinecones that are involved in this factor, we write: 

H(PC)H(PC) „ H(Pg)H(P$) Q H(P3)H(pS) 

H(P C ) = . . v - \ , H(P N ) = \V , H(P S ) = V ~' K . 

(4 e) npc )(L (°) n pc)' h[ o) nP^ 4 e) nP| 

The division by (L ( e} nP c )(L[ o) nP c ) in the first identity comes from the fact that Pf(P>) 
and Pf (P?) have edges in common at ordinates and 1. The other divisions are justified in a 
similar way. These identities, together with ([6]), give: 

H{A N )H{A S ) H(P)H{P C ) (e) „ (o) 8 (L[ o) nP»)H( P C) {L^ n Pj)H(P°) 
H(A)H(A C ) H(P»)H(PS) { 1 \ ){ \ ! (L^ nP°)H{P") (L^nP^)H(P^' 

The ratio H(P<)/(Lr°' n P c ) is the product of the edges found at non-positive ordinates in the 
horizontal matching of P< . Similarly, the ratio H(P<) / n P< ) is the product of the edges 
found at non-positive ordinates in the horizontal matching of P< . But P< C P< (see (0]) and 
its accompanying Figure [15]) , so the quotient of the two ratios is H~(P C \ P N ), the product of 
the edges found at non-positive ordinates in the horizontal matching of P< \P^. The remaining 
quotient involving P> is, symmetrically, equal to H + (P C \ P S ). This yields the result stated in 
the theorem. □ 

The refined condensation theorem specializes nicely to interleaved pinecones. 

Definition 14 (Interleaved pinecones). A closed pinecone is interleaved if, between two black 
squares in row r, one finds a black square in row r + 1 and a black square in row r — 1. 

This implies that, between two consecutive black squares in row r, there is exactly one black 
square in row r + 1, and one in row r — 1. For instance, the pinecone to the right of Figure fT71 
is interleaved. Going back to Theorem it is easy to see that for an interleaved pinecone, the 
graphs Pf \ P^ and P^ \ P^ are empty, so that p-(P c \ P N ) = H+{P C \P S ) = 1. 

Corollary 15 (The matching polynomial of interleaved pinecones). The matching poly- 
nomial of an interleaved pinecone P is related to the matching polynomial of its closed sub- 
pinecones by 



M{P)M{P C ) = ( ]J a ) M(P W )M{P E )+ I Yl 

\aG_Ei / \a£E 2 



a M{P N )M{P b ) 



where the sets E\ and are described in Theorem Moreover, the five sub-pinecones of P 
are also interleaved. 

The last statement follows from the fact that each of the five sub-pinecones can be defined as 
the largest closed pinecone containing two prescribed vertical edges. 

4.2. Special edges 

We will now simplify further the expression of Corollary [HI by assigning weight 1 to certain 
horizontal edges, called ordinary. If P is interleaved, the set of black squares of P w is obtained 
by deleting the rightmost black square in each row of P. Consequently, the rows that disappear 
when constructing P w from P are those that contain only one black square. These are the rows 
that contain a left edge contributing to the set E\. Moreover, the top and bottom rows of P 
contain exactly one black square, otherwise P would not be interleaved. Hence E\ has cardinality 
at least 2. We are going to assign weight 1 to all the edges of E\ that lie neither on the top 
segment of P nor on its bottom segment. Similarly, we will assign weight 1 to the edges of E\ 
and E\, so that the product of the edge- weights in E2 will reduce to ew. As we want to apply 
the condensation theorem iteratively, this forces us to set to 1 the weights of other horizontal 
edges, occurring for instance in the sets E\ and E\ associated to the five sub-pinecones of P. 
Iterating this procedure, we arrive at the following definition of ordinary horizontal edges — 
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those that will have weight 1. This definition is illustrated in Figure fTTl Note that it does not 
assume that the pinecone is interleaved. 

Definition 16. An even horizontal edge a, lying at ordinate r in a pinecone (that is, between 
rows r — 1 and r), is ordinary if the closest black square found in rows r — 1 and r weakly to the 
right of a is in row r — 1. Otherwise, a is said to be special. In particular, if an even edge a lies 
in the bottom segment of P, it is special. 

Symmetrically, an odd horizontal edge a, lying at ordinate r, is ordinary if the closest black 
square found in rows r — 1 and r weakly to the right of a is in row r. Otherwise, a is said to be 
special. In particular, if an odd edge a lies in the top segment of P, it is special. 




Figure 17. The ordinary edges of a pinecone. The even ones are in black, the 
odd ones in white. The pinecone to the right is interleaved. 

It is easy to check that in an interleaved pinecone, the edges of E\ and -E| are ordinary. The 
following lemma tells which edges of E\ are special. 

Lemma 17. Let P be an interleaved closed pinecone. There are exactly two left edges of P that 
do not belong to P w and are special. One of them is even, and is the lowest left edge of P. The 
other is odd, and is the highest left edge of P. 

Proof. As noted at the beginning of this subsection, the left edges of P that do not belong to 
P w are those that belong to rows containing exactly one black square. Take an even edge of this 
type. It belongs to the bottom portion of P. Figure [18] shows that it is always ordinary, unless 
it lies on the bottommost horizontal segment of P. The proof is similar for odd left edges. □ 



Figure 18. The even ordinary left edges of P \ P w . 

Lemma 18. Let P be a closed pinecone, and Q one of the five sub-pinecones P c , P w , P E , 
P N , P s . The ordinary edges of Q are exactly the ordinary edges of P belonging to Q. 

Proof. Let a be an even ordinary edge of P, lying at ordinate r. Let c be the first black square 
found in row r — 1 weakly to the right of a. By definition of ordinary edges, there is no black 
square in row r between a and c. Assume a belongs to Q and is not ordinary in Q. Since we do 
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not add squares when going from P to Q, this means that c does not belong to Q. Then there is 
no black square in row r — 1 to the right of a in Q. However, since a belongs to Q, there must be 
a black square d to the right of a in row r of Q. This square d is also in P, and to the right of c. 
But Q is defined as the largest closed subpinecone of P having certain prescribed rightmost and 
leftmost edges, so that if it contains a and c', it has to contain c as well. We have thus reached 
a contradiction, and a is ordinary in Q. 



Conversely, assume a is special in P, but ordinary in Q. The latter property implies that 
there is a black square c in row r — 1 of Q to the right of a. Of course, c also belongs to P. 
Since a is special in P, there is a black square d in row r of P lying between a and c. As Q is 
the largest pinecone containing two prescribed edges, and contains a and c, the square d must 
be in Q as well, contradicting the assumption that a is ordinary in Q. 



L I I I I I I 



in 

Q 



Of course, the proof is completely similar for odd special edges. □ 

4.3. The partial matching polynomial 

For any pinecone P, define its partial matching polynomial M(P) to be the value of M(P) 
when the weights of all ordinary edges are set to 1. We emphasize that this polynomial counts 
perfect matchings (all vertices of P belong to an edge in the matching) , but some of the edges 
have weight 1. Assume P is interleaved, and apply Corollary[TH As observed after Definition [TBI 
all the edges of E\ and E\ are ordinary, so that they have weight 1. This means that the second 
monomial occurring in the condensation formula is simply ew. Moreover, the special edges of 
Ei are described in Lemma fTTl This, combined with Lemma fl8l implies the following corollary. 

Corollary 19. The partial matching polynomial of an interleaved closed pinecone P is related 
to the partial matching polynomials of its sub-pinecones by 

M(P)M(P C ) = aa'M(P w )M{P E ) + ewM(P N )M{P s ), 

where a and a' are the highest and lowest left edges of P. 

Since the Gale- Robinson pinecones constructed in Section [3] are interleaved, we have obtained 
a combinatorial interpretation of the Gale- Robinson polynomials. 

Theorem 20. Let P{n) = P(n;i,j,k,£) be the sequence of pinecones associated with the pa- 
rameters (i,j,k,t). Let q(n) = (n;u,v) be the polynomial in u and v that counts the perfect 
matchings of P{n) according to the number of vertical edges (the exponent of v) and horizontal 
special edges (the exponent of u). Then q(n) = 1 for n < m and for n > m, 

q(n)q(n — m) = u 2 q(n — i)q(n — j) + v 2 q(n — k)q(n — £). 

This proves Theorem HU as the recurrence shows that q(n; u, v) — p(n; u 2 , v 2 ). 
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5. Perspectives 
5.1. Variations and extensions 

There is a good deal of overlap between this article and the paper by David Speyer on the 
general octahedron recurrence, of which the Gale- Robinson recurrence is a very special case [25] . 
Speyer's method allows him to construct, for each k,£) with i + j = k + £, a sequence of 
graphs having the same number of perfect matchings as the pinecones we construct. We believe 
that our graphs are the same as the ones that are given by Speyer's procedure, but we have not 
proved that this holds in general. 

One undesirable feature of our description of Gale-Robinson pinecones is that it breaks some 
of the symmetries between the parameters i, j, k, and I. Clearly, exchanging k and I reflects the 
pinecone across a horizontal line. But the convention j = min{i, j, k,£} implies that i and j do 
not play symmetric roles, nor the pairs {i,j} and {k,f\. This explains why the description of 
the bivariate polynomials q(n;u,v) of Theorem [20] is not symmetric in u and v. Perhaps some 
of this asymmetry is unavoidable, but it would be good to find a more symmetrical definition or 
else achieve some insight into why the asymmetry cannot be avoided. 

Indeed, part of the point of view that led to both this article and Speyer's is that the truly 
fundamental objects of study are functions that map a three-dimensional lattice to some ring 
and that obey the octahedron relation 

/(x + i)/(x - i) + /(x + j)/(x - j) + /(x + k)/(x - k) = 

(where x is an arbitrary vector in the lattice and i, j, k are fixed generators of the lattice) 
and more general versions of the relation that include coefficients of various kinds. There is 
no intrinsic "arrow of time" here (as there is when one thinks of running a recurrence relation 
forward from some set of initial conditions) , but some sets of initial conditions are sufficiently 
large that they allow one to reconstruct the entirety of /, and some of these subsets of the lattice 
can be viewed as "space-like", so that one can think of the reconstruction of successive slices of 
the lattice as a kind of propagation. In the fully symmetrical version, there is no reason to 
privilege one direction over its reverse, or one axis over another. 

In contrast, when one descends from this level to the more concrete world of graphs and perfect 
matchings, the symmetry appears to be broken. A full theory of the octahedron recurrence would 
incorporate graph-theoretic analogues of all the symmetries of the three-dimensional lattice; such 
an understanding is currently lacking. Just as Ehrhart theory for enumeration of lattice-points in 
polytopes can best be understood in a context that includes inside-out polytopes [2], the theory 
of Aztec diamonds, crosses-and-wrenches, and pinecones requires notions of geometric graphs in 
which combinatorial parameters that are ordinarily required to be positive can take on negative 
values as well. (E.g., one needs a theory in which the notion of an Aztec diamond of order 4 
and an Aztec diamond of order —4 enter on an equal footing, and the latter graph turns out to 
be essentially the same things as an Aztec diamond graph of order 3.) As a hint of what such a 
theory might look like, the interested reader should look at [19] and pQ. 

The bivariate polynomials p(n; w; z) studied in Section 0] generalize Gale- Robinson numbers. 
A different extension of these numbers comes from replacing the initial conditions (a string of 
m l's) by generic initial conditions (indeterminates xq through x m -i). Here again, Fomin and 
Zelevinsky proved algebraically, and Speyer proved combinatorially, that the rational functions 
one obtains are Laurent polynomials in xq, ■ ■ ■ ,x m -\. Speyer's work shows that these variables, 
in contrast to the formal coefficients w and z mentioned above, are most naturally viewed 
as being associated with the faces of a graph, rather than its edges. So there should be a 
way to associate these m variables with the faces of our pinecones and use them to assign 
weights to the perfect matchings so that the weight of each perfect matching of a pinecone is a 
Laurent monomial in Xo, . . . , x m -i. Indeed, there should be an extension of Theorem [20l that 
describes the Laurent polynomials that arise from setting a(n) — x n for < n < m — 1 and 
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a(n) = (wa(n — i)a(n — j) + za(n — k)a(n — £))/a(n — m) for n > m, and in particular identifies 
each Laurent monomial in a(n) as the weight of a perfect matching of P(n; k, £). 

Most of the work of this article was done in 2005 and 2006, as the study of cluster algebras 
was beginning its (still continuing) outward explosion, so there are now other approaches to 
proving positivity results that have some overlap with the approach taken here. In particular, it 
is possible that pinecones graphs can also be viewed as Aztec diamond graphs with defects, in 
the manner of [9]. 



5.2. Random pinecone matchings 

A rather different direction that might be studied is the "typical" behavior of perfect match- 
ings of large pinecones. Figure [19] shows two tilings associated with matchings of Somos-4 
pinecones. (Here we make use of the standard duality between a tiling of a polyomino by 
dominos and a perfect matching of the dual graph of the polyomino, in which vertices cor- 
respond to cells of the polyomino and edges correspond to pairs of adjoining cells, i.e. legal 
positions of a domino in a tiling.) The first one corresponds to n = 26 (that is, to a perfect 
matching of the graph P(26; 3, 1, 2, 2)), the second one to n = 50. Both were chosen uni- 
formly at random from the set of all perfect matchings of that graph. These examples were 
produced using Propp and Wilson's papers on "exact sampling" [HI [2Q] which show how the 
method of "coupling from the past" permits one to generate random perfect matchings of bipar- 
tite planar graphs. Indeed, this algorithm was incorporated into a program called vaxrandom 
that accepts a VAX-file as input and produces a perfect matching of the associated graph as 
output, or rather, the dual picture of a domino tiling of a region. The source code for the 
program is contained in the files |http : //jamespropp . org/tiling/sources/vaxrandom. c| and 



http://jamespropp.Org/tiling/sources/allocate.h, and information on the program's use 
can be found at |http : // j amespropp . org/ 1 iling/ doc/ vaxrandom . html ) 
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The reader will quickly notice that in both of these random tilings, the randomness is not 
spatially distributed in a uniform manner. Near the boundary, there is a good deal of order, with 
tiles lined up the same way as their neighbors; only in the interior does one find random-looking 
behavior. 

This phenomenon is not specific to pinecones, but has been observed for a wide variety of 
two-dimensional tiling models over the past decade, from [7] and [5] to [16]. The most-studied 
case is the Aztec diamond graph ((i,j,k,£) = (1,1,1,1), in our notation); in this case, it has 
been shown that in a suitable asymptotic sense there is a sharp boundary between the part of 
the tiling that is random and the part that is orderly, and that this boundary is (asymptotically) 
a perfect circle. A similar sort of domain-boundary is visible in Figure [191 assuming that the 
theory for pinecones is analogous to the theory for Aztec diamond graphs, it would be interesting 
to know the asymptotic shape of the domain-boundary for (i, j, k, ^)-pinecones asrn oo. 

One interesting feature of Gale-Robinson pinecones is that we can write the definition in a 
way that makes sense even when the parameters (i, j, k,£) cease to be integers. Formula |(2]) can 
be rewritten as 

U(t, r, c) = 2c + r - 3 - 2 [fie + nr + i - t\ , 

(7) 

L(t, r,c) = 2c + r - 3 - 2 |_£tc + Ar + i - t\ . 

where l = i/j, k = k/j, A = £/j, u — m/j = l+1 = k + X, and t = (n + So there is 

a sense in which all the pinecones discussed in this article are part of a four parameter family, 
parametrized by i, k, A and t. Of course, the graphs do not vary continuously in these variables 
(being discrete elements in a countable set, namely the set of all finite graphs, how could they?), 
but this parametrization seems likely to be natural for some purposes, e.g., the study of random 
perfect matchings of pinecones. (It is to be expected that a coherent limit-law with t — > oo 
will prevail for any fixed choice of (l,k, A), whether or not t, k, and A are rational.) It should 
be noted, incidentally, that if one chooses parameters (i,j,k,i) with a greater common divisor 
d > 1, the sequence of pinecones one gets from our construction is the same as the sequence of 
pinecones that one gets from the parameters (i/d,j/d,k/d,£/d), except that each pinecone in 
the latter sequence is repeated d times in the former sequence; this observation follows easily 
from the i, k, A formulation of the definitions of U(-) and L(-). 

5.3. Closed-form expressions 

One feature common to sequences satisfying three-term or four-term Gale-Robinson recur- 
rences is that the terms grow at quadratic exponential rate. Indeed, it is easy to verify, from the 
discussion of pinecones, that in the infinite sequence of graphs associated with any particular 
three-term Gale-Robinson recurrence, the nth graph has 0(n 2 ) vertices, with each vertex having 
degree at most 4. It follows from this that Gale-Robinson sequences have at most exponential- 
quadratic growth; that is, the nth term is bounded above by C" for all sufficiently large C . 
In some cases, an exact formula is possible; we have already mentioned the 'Aztec diamond 
case" i = j = k = £— 1, and in the case (i, j,k,£,m) = (6, 1,4,3,7) there is an exact formula 
for a(n) of the form 2 e2 (")3 e3< -™- ) where the exponents e2(n) and e^{n) are given by quadratic 
polynomials in n whose coefficients are periodic functions of n (we thank Michael Somos for 
bringing this special case of the Gale-Robinson recurrence to our attention, and we raise the 
question of whether there are other instances of Gale- Robinson sequences being given by simple 
exact formulas). However, in general such algebraic formulas do not exist. Instead, one must be 
content with formulas that express the nth term in terms of Jacobi theta functions. This link 
with the analytic world is what motivated Michael Somos to introduce the Somos-fc sequences 
to begin with. E.g., back in 1993, Somos announced (without proof) that the nth term of the 
Somos-6 sequence is given by/(n — 2.5,n — 2.5) where 

oo oo 

f(x,y)= Cl 4 y J2 (- 1 )" 2 E 9(ki,k 2 ,x,y), 

fc 2 — — oo ki — — oo 
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g(k 1 ,k 2 ,x,y) = c^c* 2 ^ 1 ** cos(c 6 kix + c 7 k 2 y), 



Cl = 


0.875782749065950194217251... 


C2 = 


1 .084125925473763343779968. . . 


C3 = 


0. 1 14986002186402203509006. . . 


c 4 = 


0.077115634258697284328024... 


C5 = 


1 .180397390176742642553759. . . 


C6 = 


1 .508030831265086447098989. .. 


C7 = 


2.551548771413081602906643.. 



(See http://jamespropp.org/somos/ellipticfor a similar but simpler formula for the Somos- 
4 sequence.) However, as far as we are aware, nobody has proposed (or even conjectured) a 
fully general analytic formula for the terms of sequences satisfying three-term Gale-Robinson 
recurrences. A more detailed discussion of the analytic properties of such sequences can be 
found in [15], which also gives some of the history of these sequences. 

It is worth mentioning that for the Somos-4 sequence, there exists a unique constant c such 
that s(n) (the nth term of the sequence) is on the order of c™ , but that the behavior of s(n) / c n 
is oscillatory; see |http : / /j amespropp . org/ somos /elliptiq Let us also mention a recent paper 
by Xin [57] where the Somos-4 numbers are expressed as determinants of Hankel matrices with 
integer coefficients. 

5.4. Analogy with the KP hierarchy 

We conclude with some remarks (based on some unpublished remarks of Andrew Hone) about 
the analogy between Somos sequences and the like and the hierarchy of solutions to an integrable 
PDE like the KdV equation, followed by our own speculation about a direction for further study 
that the analogy might suggest. 

The equation 

u X xx + 6uu x +u t = 0, 

where u = u(x,t) is the function we want to solve for and subscripts indicate partial differen- 
tiation (e.g., u xxx = is known as the KdV equation, and has played a crucial role in the 
modern theory of partial differential equation, as part of a large family of equations with related 
properties (the "KP hierarchy"). If one sets u — 2(d x ) 2 \ogF one can rewrite the PDE in the 
compact form 

(D X D X D X D X + D x D t )(F ® F) = 
where D x and D t are the "Hirota D-operators" acting on tensor-pairs of functions via 

D x (f(x,t) ®g(x,t)) = (d Xl - d X2 )f(x 1 ,t)g(x 2 ,t)\ Xl=X2=x 

and 

D t (f{x,t) ®g(x,t)) = (d tl - d t2 )f{x,t 1 )g(x,t 2 )\t 1 =t 2 =t- 
(Note that in the literature on KdV, this tensor product is traditionally written as / • g rather 
than f®g and is called the "dot-product", but it is a tensor product, not an inner product). More 
generally, the bilinear method is the trick of rewriting PDEs in the form P(D X , D y , . . .)(F®F) — 
0. Hirota operators are antisymmetric, so we can think of them as actions on the antisymmetric 
square of a vector space of functions. For more on the Hirota method, see e.g. |14| . 

Analogously, if we take V to be the vector space of real- (or complex-) valued bilaterally 
infinite sequences (. . . ,sq, . . .), we may define, for every pairs of integers a bilinear shift 
operator V ® V — > V sending (s n )-oo ® (in)^oo to (sn+itn+j)™^ (the sequence whose nth term 
is s n+i t n+ j for all n 6 Z). These operators, graded by i+j, generate a graded ring of bilinear 
shift-operators, and the Somos sequences and Gale-Robinson sequences are special instances of 
sequences (sn)^^ for which the tensor-square (sn)^oo ® ( s n,)^oo li es m the kernel of a particular 
bilinear operator. It has been noticed that for a typical Somos or Gale-Robinson sequence, the 
tensor-square of the sequence, in addition to being annihilated by the "defining" bilinear operator, 
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is annihilated by infinitely many others as well. In fact, there is more than just an analogy at 
work here: each GR recurrence can be written in terms of Hirota differential operators by taking 
exponentials (see [23]). Hopefully, by combining algebraic, analytic, and combinatorial tools, 
future researchers will shed some light on this intriguing phenomenon. 
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